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Abstract 

We show that suppression of the baryon energy gaps, caused by the relative motion of superfluid 
and normal liquid components, can substantially influence dynamical properties and evolution of 
neutron stars. This effect has been previously ignored in the neutron-star literature. 
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I. INTRODUCTION 

According to numerous microscopic calculations (e.g., [l, 2] and references therein), nu- 
cleons and hyperons in the internal layers of neutron stars (NSs) can become superffuid at 
temperatures T < 10 s 4- 10 10 K. Superfluidity has a strong impact on the thermal evolution 
of NSs, their oscillations, and (most probably) leads to such observational phenomena as 

n n n n 

glitches [3j and pulsar spin precession [4], |5|. Recent real-time observations [6| of a cooling 
NS in Cassiopea A supernova remnant give strong arguments that the star has superfluid 
core [7, 8]. 

The aim of this short note is to point out the importance of one effect related to super- 
fluidity of baryons in NSs that has usually been ignored in the NS literature. In Sec. II we 
outline the effect. In Sec. Ill we demonstrate its efficiency In Sec. IV we discuss possible 
consequences for the physics of NSs and in Sec. V we conclude. We use the system of units 
in which fcs = h — 1. 



II. A SIMPLE PROBLEM AND THE PROPOSED EFFECT 



Let us consider a degenerate Fermi-liquid composed of identical particles. Assume that 



they interact through a weakly attractive potential so that BCS theory [9] is applicable. 
Assume also that they pair (become superfluid) in the spin-singlet l S state at temperatures 
T below some critical temperature T c . The role o 
Fermi-liquid is played by Bogoliubov excitations 
written in a reference frame in which the mean (hydrodynamic) velocity V q of Bogoliubov 
excitations vanishes, V q = (i.e., normal liquid component is at rest). 

In the absence of superfluid current (when the superfluid velocity V s = 0) the energy E p 
of a Bogoliubov excitation with momentum p near the Fermi surface can be written as 



' elementary excitations in such superfluid 
lo| . In what follows, all equations will be 



E P = ^vU\p\-Pf) 2 + A*, (1) 

where v-p and pp are the Fermi-velocity and Fermi-momentum, respectively; and A is the 
energy gap, given by the standard equation 

l = -V.E^. (2) 
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where Vq is the (constant) pairing potential and 

f P = e E P /T +1 ( 3 ) 

is the Fermi-Dirac distribution function for Bogoliubov excitations. 

If, however, the superfluid current is present (V s 7^ 0) then fermions pair with momenta 
(— p + Q,p + Q) rather than with (— p,p), and the total momentum of a Cooper pair is 

2Q = 2mV s . (4) 

What will be the equation for the gap? The answer can be found in the paper by Bardeen 



ll| and is well known in the physics of superconductors. Now, instead of Eq. (j2J), one should 



write 

1 = _ Vq J2 1 _^±9_I^±Q. (5 ) 



2E P 

p v 



Here F p+ q is the distribution function for Bogoliubov excitations with momentum (p + Q) 
in the system with non-zero V s , 

= e W r + l' (6) 



where 



« — + E P (7) 

m 



is the energy of a Bogoliubov excitation with momentum 
Q <C Pf which is true in all interesting cases (see, e.g., [12 



p — Q ) . In Eq. ((7j) we assumed 



Eq. (jSJ) can be written in terms of the quantity Ao, which is the energy gap at T = and 
Q = 0. It satisfies Eq. (j2]) with / p = 0. Using it, one can present Eq. (JS]) in the form 

The solution to this equation gives the gap A as a function of T and Q — \Q\. 
First consider two limiting cases in which A(T, Q) vanishes, 
(z) if Q = then A = at 

T = T C ^ 0.567A (the well known BCS result); (9) 

fill if T = then A = at 
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FIG. 1: (color online) The energy gap A (in units of Aq) versus Q = mV s [in units of Q cr o, see Eq. 
(033)] for a set of temperatures T/T c = 0.1, 0.4, 0.6, 0.75, 0.85, and 0.95. 

The latter result is less known but can be found, e.g., in [15]. Notice, that the well-known 
Landau criterion for superfluidity breaking gives Q^Q Ildau ' ) = A m/p F and is not accurate 
for a superfluid Fermi-liquid. 

Some numerical solutions to Eq. (jSJ) are presented in Figs. 1 and 2. Fig. 1 shows the gap 
A(T, Q) [in units of Ao] versus momentum Q [in units of Qcro] for a set of temperatures 
T/T c = 0.1, 0.4, 0.6, 0.75, 0.85, and 0.95. One sees that A is quite sensitive to variation 
of Q = mV s as long as T > 0.1T C . Another important conclusion that can be drawn from 
Fig. 1 is that (for a given T) the maximum critical momentum Q CT strongly depends on 
temperature. 

Fig. 2 illustrates this point more clearly. In the left panel we plot Q CT (in units of Q ct q) 
versus T (in units of T c ). The right panel shows the same dependence Q cr (T) but with Q cr 
measured in units of 

e'r'p-) = § (id 

We see that Q cr changes with T in such a way that Q cr (T) /Q& (T) is roughly constant. 

Therefore, the energy gap A can be a strong function of the momentum Q = mV s or, in 
an arbitrary frame, a strong function of the difference m (V„ — V q ) = in AV. We will refer 
to this effect as to the 'AV-effect'. The critical value AV CT (T) of AV = |V S — V q |, at which 
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FIG. 2: (color online) Left panel: Q cr (in units of Q C ro) versus T (in units of T c ). Right panel: 



The same as in the left panel but Q cr is in units of Q, 



(app) 



superfluidity dies out, is easily estimated by taking Q cv ~ Qcr PP ' ) - Then, from Eq. (ITT]) , we 
obtain 



AV CT (T) ~ 10 7 



A(T, 0) 
10 9 K 



^o\ 1/3 _i 
— I cm s , 
n / 



(12) 



where A(T, 0) is measured in Kelvins; n = 0.16 fm 3 is the nucleon density in atomic 
nuclei; n = pp/(3ii 2 ) is the particle number density. 



III. IMPORTANCE OF THE AV-EFFECT FOR NEUTRON STARS 

If the difference AV between the baryon superfluid velocities and a normal velocity is 
comparable to AV cr , then the baryon energy gaps can be substantially reduced. A few in- 
teresting consequences of this 'dynamical reduction' of the gaps are discussed in the next 
section. Here we illustrate possible importance of the AV-effect by considering radial os- 
cillations of a nonrotating superfluid NS whose core is composed of neutrons, protons, and 
electrons. The main question is at what oscillation amplitude AV^ becomes comparable to 
AV cr ? 

For simplicity we (i) assume that neutrons pair in the spin-singlet ( 1 5'o) state [rather 
than in the triplet ( 3 -P2) state] and (ii) neglect the Landau quasiparticle interaction between 
quasinucleons when calculating A n (T, V sn — V q ) [here and below the subscripts n and p refer 
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FIG. 3: (color online) Left panel: Nucleon critical temperatures T cn and T cp versus density p. 
Right panel: Redshifted critical temperatures and T™ versus radial coordinate r (see text for 
details). 



to neutrons and protons, respectively]. In [16| we outline the calculation of A n and A p as 



functions of (V sn — V q ) and (V sp — V q ) with allowance for interactions between quasiparticles. 
The NS model used here and all the microphysics input are essentially the same as in 



171 ]; we refer the reader to that work for more details. In particular, we consider the star 
of gravitational mass M = 1.4M , circumferential radius R = 12.2 km, central density 
p c = 9.26 x 10 14 g cm -3 , and adopt the APR EOS in the NS core |18]. The model of nucleon 



rr 

[17] • 



superfluidity employed here coincides with the model 3 of Ref. [17] and is shown in Fig. 3. 

The left panel of Fig. 3 presents nucleon critical temperatures T cn and T cp versus density 
p in the NS core, the right panel demonstrates the red-shifted critical temperatures = 
T cn e V//2 and = T cp e u ^ 2 (y is the metric function) versus radial stellar coordinate r (in 
units of R). The redshifted proton critical temperature is taken to be constant = 
5 x 10 9 K; the redshifted neutron critical temperature varies with r and has maximum 
T™ mai = 6 x 10 8 K in the stellar centre. In the right panel of Fig. 3 we hatch the region 
occupied by the neutron superfluidity at a redshifted stellar temperature T°° = Te^ 2 = 
4 x 10 s K. 

To model oscillations of superfluid NSs one has to use the hydrodynamics of mixtures 
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FIG. 4: (color online) Amplitudes of the eigenfunctions AV n (solid lines) and the critical velocities 
AV cr (dashes) versus r/R for the four temperatures T°° = 3.0 x 10 7 K (black lines), 8.0 x 10 7 K 
(red lines), 2.0 x 10 8 K (blue lines), and 5.0 x 10 8 K (violet lines). To plot AV n we assumed that 
the energy of oscillations is E mec ^ = 10 47 erg. The vertical dotted lines show the (temperature- 
dependent) boundaries between the inner superfluid and the outer normal regions. See text for 
details. 



of superfluid Fermi-liquids 



19H22j. The imp ortant parameter of such hydrodynamics is 



the so called entrainment matrix 
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191 (or relativistic entrainment matrix K 



141]). which is very temperature-dependent [12|, PJj. As a consequence, the eigenfrequencies 



ik 



and eigenfunctions of oscillating superfluid NS are also temperature-dependent |17l. |21|. |23|. 
Below we consider the first radial oscillation mode of a superfluid NS (see [17|, particularly 
figure 3 there). 

Figure 4 shows the amplitude of the eigenfunction AV n = |V sn — V q | and the critical veloc- 
ity AV CI as functions of r (solid and dashed lines, respectively; both in units of 10 7 cm s _1 ). 
We plot AV n and AV cr for four redshifted stellar temperatures: T°° = 3.0 x 10 7 K (black 
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FIG. 5: (color online) Neutron energy gap A(T, AV n ) (in units of 10 9 K) versus r/R for two 
temperatures T°° = 8 x 10 7 K (left panel) and T°° = 2 x 10 8 K (right panel) and some oscillation 
energies E mec h (indicated in the figure) . Vertical dotted lines show r at which neutron superfluidity 
disappears (A n = 0). The larger -Emech the smaller the superfluid region. See text for details. 

lines), 8.0 x 10 7 K (red lines), 2.0 x 10 8 K (blue lines), and 5.0 x 10 8 K (violet lines). The os- 
cillation frequencies u of the first radial mode for such temperatures are u/ (10 4 s -1 ) ~ 1.702, 
1.702, 1.064, and 0.516, respectively. 

The vertical dotted lines in Fig. 4 indicate (temperature-dependent) boundaries between 
the neutron superfluid region and the outer normal region with nonsuperfluid neutrons. In 
the normal region the functions AV n and AV CI are not defined. The oscillation energy of the 
star is -E mec h = 10 47 erg. For a nonsuperfluid NS this energy corresponds to an oscillation 
amplitude 

e = lim ^ « 4.4 x 10" 4 , (13) 

r->0 r 



where £(r) is the Lagrangian displacement 24 1. 



It follows from Fig. 4 that AV n can substantially exceed the critical values AV CI , so that 
superfluidity is destroyed by oscillations in the large part of the stellar core (see the violet 
and blue curves). This means that the AV-effect can greatly influence (or even drive) the 
dynamics of NSs already at rather modest oscillation amplitudes. 

This point is additionally illustrated in Fig. 5, where we plot the neutron energy gap 



8 



TABLE I: Mechanical energy E mec ^ an d the corresponding amplitude of oscillations e, defined by 
Eq. PJ. 



^mech, erg 


e 


0.0 


0.0 




10 46 


1.4 x 10" 


-4 


5.0 x 10 46 


3.1 x 10" 


-4 


10 47 


4.4 x 10" 


-4 


5.0 x 10 47 


9.7 x 10" 


-4 


10 48 


1.4 x 10" 


-3 


5.0 x 10 48 


3.1 x 10" 


-3 



A n (T, AV n ) versus r/R for two temperatures, T°° = 8.0 x 10 7 K (left panel) and T°° = 
2.0 x 10 8 K (right panel), and a set of oscillation energies -E mec h- In the left panel A n (T, AV n ) 
is shown for E mcch = , 5.0 x 10 47 , 10 48 , and 5 x 10 48 erg; in the right panel A n (T, AV„) is 
shown for -E mcc h = 0, 10 46 , 10 47 , and 10 48 erg. The oscillation amplitudes e [given by Eq. 
([13"]) ] for these oscillation energies are presented in Table. 

Notice that in each panel of Fig. 5 the curves are plotted using the eigenfunctions AV n (r), 
which differ from one another only by normalization (by the value of E mec h). For -E mec h — 
10 erg these eigenfunctions have already been presented in Fig. 4 (see the red and blue solid 
lines; the red line corresponds to T°° = 8.0 x 10 7 K, the blue line - to T°° = 2.0 x 10 s K). 

If -E m ech = (no oscillations; see the solid lines in both panels of Fig. 5) the gap A n is 
unaffected by AV n and is entirely determined by the dependence of T cn on r (see Fig. 3). 
The vertical dotted lines in Fig. 5 indicate boundaries between the inner superfluid and the 
outer normal regions; these boundaries depend on -E mech . Obviously, the higher -E mec h, the 
larger AV n (r), and, correspondingly, the smaller the superfluid region and A n . One sees 
that the gaps are very sensitive to variation of AV n . 
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IV. POSSIBLE APPLICATIONS 



As follows from the consideration of the previous section, the AV-effect can operate at 
not too small oscillation amplitudes. All interesting consequences of this effect are related 
to the reduction of baryon gaps. Let us list some of them: 

(1) The reduction of the gaps influences the entrainment matrix p^ 12], which depends 
on them. As a result, p^ will become a non-linear function of the oscillation amplitude. This 
will (i) make the oscillation equations nonlinear and hence (ii) affect the eigenfrequencies and 
eigenfunctions of oscillating NS. Moreover, this will (Hi) influence the dissipation processes, 
because bulk viscosity terms explicitly depend on p^. In a rotating star the change (decrease) 
of the element Pnp of the entrainment matrix will, in addition, (iv ) decrease the mutual 



friction force 



25] . which is proportional to p np . (We emphasize that the dependence of 



p np on T and on AV n and AV P is a very important effect for mutual friction and related 
phenomena, which has been neglected in the literature.) 

How to calculate the entrainment matrix p^ taking into account the AV-effect? A direct 
calculation is difficult (but one can make it in a manner similar to how it was done in 12]). 
A good approximation for p^ could be to take the result of 12] and calculate p^ with the 
velocity-dependent gaps from Sec. II instead of the gaps A n (T, 0) and A P (T, 0). 

(2) Another important consequence of the AV-effect is its impact on kinetic coefficients 
of NS matter, in particular, on the bulk and shear viscosities. 

(i)Bulk viscosity. There are four bulk viscosity coefficients in the npe- matter of NSs 



221 ] . All of them are generated by nonequilibrium beta-processes (direct or modified Urea 



reactions) and depend on the difference Ar between the direct and inverse reaction rates. 

unction of T, A n , A p , and of the imbalance of chemical 
271 ]. where lu is the chemical potential for particle species 
, that if 5p > max{A n , A p } then, even if 



Ar is generally a complicated 
potentials 5p = p n — p p — p 
i — n, p, e. Recently it has been shown 



26 



T <C A n and/or A p , the bulk viscosity is not suppressed by the nucleon superfluidity and 
can be very efficient. It seems that the AV-effect of the reduction of the energy gaps A n 
and A p by relative motion of superfluid and normal component is complementary to the 



effect considered in 



. Both effects act in unison to increase the bulk viscosity coefficients, 



and they are of comparable strength. Notice, however, that the effect of Ref. 28] can only 
affect the bulk viscosity coefficients, while the applicability range of the AV-effect is wider; 
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it directly influences the baryon energy gaps and thus all dynamics of NSs. 

(ii) Shear viscosity. Neglecting entrainment between baryon species (p np = 0), the shear 
viscosity 77 can be calculated in the same fashion as was done, e.g., in [29] (the results 
will be the same). The only difference is that one should use the velocity-dependent gaps 
Aj(T, AVi) instead of Aj(T, 0) [i = n, p\. It is interesting that the AV-effect can both 
increase or decrease the shear viscosity. For example, the electron shear viscosity r\ e decreases 
with increasing AV P (that is, with reducing A p , because the electrons are better screened 
by protons when A p is large 29]). On the other hand, the neutron shear viscosity r\ n can 
either decrease or increase with growing AV n and AV p . The behaviour of rj n in that case 
is determined by the competition of two effects: by the increase of the normal density of 
neutron Bogoliubov excitations p qn and by the reduction of the neutron mean free path A 
due to more frequent collisions with neutron and proton Bogoliubov excitations (note that 
rj n can be estimated as r\ n ~ p qn v-p n A, where vp n is the neutron Fermi momentum). Similar 
effects were carefully analyzed in 30[ in application to the neutron thermal conductivity. 

An entrainment between neutrons and protons will strongly modify the derivation of 
the neutron shear viscosity, even neglecting the AV-effect. The main difference will be 
the equilibrium Fermi-Dirac distribution function for neutron Bogoliubov excitations in a 
system with superffuid currents. This function was first obtained in [12J [see equation (28) 
there]; it is very different from the standard expression, valid when p np = 0. To our best 
knowledge, a derivation of r] n in a system with entrainment has not been attempted in the 
literature. 

(3) Finally, there is a number of important consequences of the fact that the relative 
velocity AV between the superfluid and normal liquid components cannot be too large in a 
stationary rotating NS. Here we present two of them. 

(i) It is generally accepted that neutron vortices are pinned to atomic nuclei in the NS 
crust (or to magnetic flux tubes in the NS core). At some critical AV they can unpin from 
the nuclei (or from magnetic flux tubes). However, in some models (e.g., 31] ) pinning is so 
strong that the critical relative velocity can be as high as 10 6 -j- 10 7 cm s" 1 . These values are 
close to AVr [see E q. (fl2l) |. Thus, the AV-effect can be very important for such models. 



(ii) In Refs. m, 



a two-stream instability is discussed that can be triggered once the 



relative velocity AV reaches some critical value. According to 



32 



33|, this critical value 



is of the order of the sound speeds, i.e., it is much greater than the typical AVr, at which 
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superfluidity completely disappears [see Eq. ( TT2|) ]. In other words, it is not very probable 
that this instability is realized in NSs. 

V. CONCLUSION 

The baryon energy gaps depend on the relative velocity between the superfluid and normal 
components (AV-effect). We propose, for the first time, that this effect may have a strong 
impact on the dynamical properties of NSs. We illustrate this point by considering radial 
oscillations of an NS with superfluid nucleon core and a nonsuperfluid crust. However, we 
stress that the AV-effect should be equally important in the crust of NSs where superfluid 
neutrons are present, as well as in the interiors of hyperon and quark stars. Although we 
discussed some immediate applications in Sec. IV, it is clear that more efforts are needed to 
analyze all possible consequences of this effect on the evolution of NSs. 
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